The dual Meissner effect in SU(2) Landau gauge 
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Abstract. The dual Meissner effect is observed without monopoles in quenched SU(2) QCD 
with Landau gauge-fixing. Abelian as well as non-Abelian electric fields are squeezed. Magnetic 
displacement currents which are time-dependent Abelian magnetic fields play a role of solenoidal 
currents squeezing Abelian electric fields. Monopoles are not always necessary to the dual Meissner 
effect. The squeezing of the electric flux means the dual London equation and the massiveness of 
the Abelian electric fields as an asymptotic field. The mass generation of the Abelian electric fields 
is related to a gluon condensate < A" A" >^ of mass dimension 2. 



INTRODUCTION 

To understand color confinement mechanism is still an unsolved important problem. 
The dual Meissner effect is believed to be the mechanism[l, 2]. However what causes 
the dual Meissner effect is not clarified. A possible candidate is magnetic monopoles 
which appear after projecting SU (3) QCD to an Abelian U( 1 ) 2 theory by a partial gauge 
fixing[3]. If such monopoles condense, color confinement could be understood due to 
the dual Meissner effect. Actually an Abelian projection adopting a special gauge called 
Maximally Abelian gauge (MA)[4, 5] leads us to interesting results[6, 7] which support 
the idea of monopoles after an Abelian projection. See Fig.l in which the dual Meissner 
effect due to monopole currents as a solenoidal current is seen beautifully [8]. 




FIGURE 1. Abelian electric field flux and monopole currents in MA gauge[8]. 
Now a natural question arises. What happens in other general Abelian projections 
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FIGURE 2. Abelian and non-Abelian E electric field profiles in Landau gauge. W(R x T = 6 x 6) 
is used. 



or even without any Abelian projection? For example, consider an Abelian projection 
diagonalizing Polyakov loops. Monopoles exist in the continuum limit at a point where 
eigenvalues of Polyakov loops are degenerate[3]. But it is easy to show that such a point 
runs only in the time-like direction. This means only time-like monopoles which do 
not contribute to the string tension exist in the Abelian projection[9]. Discuss another 
simple case of Landau gauge. There vacuum configurations are so smooth and it is easy 
to check numerically that no monopoles coming from singularities exist. Without space- 
like monopoles or monopoles themselves, monopole condensation could not occur. We 
have to find another confinement mechanism. 

In this note, we show that the dual Meissner effect in an Abelian sense works good 
even when monopoles do not exist, performing Monte-Carlo simulations of quenched 
SU (2) QCD with Landau gauge fixing. Instead of monopoles, time-dependent Abelian 
magnetic fields regarded as magnetic displacement currents are squeezing Abelian elec- 
tric fields. The dual Meissner effect leads us to the dual London equation and the mass 
generation of the Abelian electric fields which suggests the existence of a dimension 
2 gluon condensate. Our present numerical results are not perfect, since the contin- 
uum limit, the infinite- volume limit and the gauge-independence are not studied yet. 
Moreover our discussions use Abelian components only on the basis of not yet clari- 
fied assumption that Abelian components are dominant in the infrared QCD (Abelian 
dominance[6, 7, 10, 11]). Nevertheless authors think the present results are very inter- 
esting to general readers, since they show for the first time the Abelian dual Meissner 
effect is working in lattice non-Abelian QCD without resorting to monopoles coming 
from a singular gauge transformation[12]. The gauge adopted here is the simplest one 
only to get smooth configurations. The present results hence suggest the Abelian dual 
Meissner effect is the real universal mechanism of color confinement which has been 
sought for many years. Moreover the relation of the Abelian dual Meissner effect with 
the dimension 2 gluon condensate sheds new light on the importance of the gluon con- 
densate [13, 14, 15, 16, 17]. Detailed studies, a direct proof of Abelian dominance and 
extensions to SU (3) and full QCD in zero and finite-temperature cases are in progress. 



SIMULATION DETAILS 



As a lattice quenched SU(2) QCD, we use an improved gluonic action found by 
Iwasaki[18]. 

S = /3 {co^Tr(plaquette) + c\^Tr(rectangular) } , (1) 

with which better scaling behaviors of physical quantities are expected. The mixing 
parameters are fixed as cq + Sc\ = 1 and c\ = —0.331. 

To measure correlations of gauge-variant electric and magnetic fields directly, we 
adopt a simplest gauge called Landau gauge which maximizes Y, s ,n Tr[Up,(s) + U^(s)]. 
After the gauge fixing, we try to measure electric and magnetic flux distributions by 
evaluating correlations of Wilson loops and field strengths. For comparison, we also 
use MA gauge where £ s , ^ Tr[L^ (5)036^ (s) 03] is maximized. To get a good signal 
to noise ratio, the APE smearing technique[19] is used when evaluating Wilson loops 
W(R,T) = W° + iW a o a : 

Ui(s) — N^Ui{s) + aJ^Uj{s)Ui{s + j)u](n + 1) J , (2) 

where is normalization factor and a is a free parameter. We have used a = 0.2 and 
N = 80. 

Measurements of the string tension make us fix the scale when we use -Jci = 440MeV. 
We adopt a coupling constant j8 = 1.2 in which the lattice distance a(/3 = 1.2) is 
0.07921 (22) [fm]. The scale is chosen only because we compare our results with those 
studied extensively in MA gauge[20, 8]. The lattice size is 32 4 and after 5000 ther- 
malization, we have prepared 5000 thermalized configurations per each 100 sweeps for 
measurements. 

Non-Abelian electric and magnetic fields are defined from lxl plaquette U ilv {s) = 
U° v + iU^ v o a as done in Ref.[21]: 

Et(s) = \{U a Ak {s-k) + U a Ak {s)) 

B a k (s) = l -e k UUf m {s-t-m) 

+U? m (s-t) + U? m (s-m) + Ut m (s)) 

We also define Abelian electric (£^.) and magnetic fields (5^.) similarly using Abelian 
plaquettes 9 c ^ v (s) defined through link variables 6£(s): 

e?vM = W + e?(*+A) - o a ^(s+v) - ez( s ) (3) 



where 6„(s) is given by U^(s) = exp(i6£(s)o a ). In MA gauge, the Abelian link variables 
(s) are defined by a phase of the diagonal part of the non- Abelian link field. 

Ufa) = ^l-Ms)\ismGf? A (s). 

Since the off-diagonal part \cfx (s) | is small[6], 0^ A (s) ~ 0^ (s) in MA gauge. As a source 
corresponding to a static quark and antiquark pair, we use here only non- Abelian Wilson 
loops. 



RESULTS 

First we show in Fig. 2 Abelian and non-Abelian electric flux profiles around a pair 
of static quark and antiquark in Landau gauge. The profiles are mainly studied on the 
perpendicular plane at the midpoint between the quark pair. Note that electric fields 
perpendicular to the QQ axis are found to be negligible. It is very interesting to see from 
Fig. 2 that Abelian electric field E& z simply defined here is squeezed also. Moreover the 
squeezing of E Az is stronger than that of non-Abelian one E z . To know how squeezing 
of the Abelian flux occurs seems hence essential. 

Let us discuss from now on flux distributions of Abelian fields alone. It is checked 
numerically that there are no DeGrand-Tous saint monopoles[22]. See Fig. 3 in which 
histograms of Abelian field strength 9^ are plotted in Landau gauge (Left) and MA 
gauge (Right). 
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FIGURE 3. Histograms of MV in Landau (Left) and MA (Right) gauges 
Hence the Abelian fields satisfy kinematically the simple Abelian Bianchi identity: 

VxE A 1 = d 4 B%, W-B a A = Q. (4) 

In the case of MA gauge, there are additional monopole current (k, ka) contributions: 

VxF = d 4 B MA + k, V-B MA = k 4 . (5) 




FIGURE 4. Magnetic displacement currents in Landau gauge as a solenoidal current. 



Here E MA and B MA are defined in terms of plaquette variables Qff A (s) (mod lit) which 

are constructed by 0^ A (s). 

The Coulombic electric field coming from the static source is written in the lowest 
perturbation theory in terms of the gradient of a scalar potential. Hence it does not 
contribute to the curl of the Abelian electric field nor to the Abelian magnetic field in the 
above Abelian Bianchi identity Eq.(4). The dual Meissner effect says that the squeezing 
of the electric flux occurs due to cancellation of the Coulombic electric fields and 
those from solenoidal magnetic currents. In the case of MA gauge, magnetic monopole 
currents k play the role of the solenoidal current[23, 20, 8]. 
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FIGURE 5. Curl of Abelian electric fields and magnetic displacement currents around a static quark 
pair in Landau (Left) and in MA (Right) gauges. Monopole currents are also plotted in MA gauge. 

Now what happens in a smooth gauge like the Landau gauge where monopoles do 
not exist? From Eq.(4), only regarded as a magnetic displacement current could 
play the role of the solenoidal current. It is very interesting to see Fig. 4 in which this 
happens actually in Landau gauge. Note that the solenoidal current has a direction 
squeezing the Coulombic electric field. Let us see also the detailed distributions shown 
in Fig. 5. The other components of the magnetic displacement current difi^r and 9^Ba z 
are not vanishing but they are much suppressed consistently with Fig. 2. In comparison, 




FIGURE 6. Abelian, monopole and photon static potentials. 



we show the case of MA gauge also in Fig. 5. Here c^-Ba^ is found to be negligible 
numerically as already expected from the works[23, 11, 20]. Instead monopole currents 
are circulating[23, 20, 8]. In this case, k r is non- vanishing, although it is also suppressed 
in comparison with k&. k z is almost zero. The authors think that non- vanishing of the 

radial and z components of d^B in Landau gauge and k in MA gauge is due to lattice 
artifacts and also due to the smallness of the Wilson loop size adopted here. It is 
interesting that the shapes of d^B^ in Landau gauge and k<h in MA gauge look similar, 
although the strengths are different. They have a peak at almost the same distance around 
0.2[fm] and almost vanish around 0.7[fm]. 



ABELIAN DOMINANCE TEST 



The reader may wonder if the above consideration of the Abelian fields alone is enough. 
Some people believe that the non-perturbative confinement problem in infrared QCD 
could be understood in terms of Abelian quantities [6, 7, 24]. 

First we evaluate the Abelian and monopole contributions to the string tension in MA 
gauge using the Iwasaki improved action (1). It is plotted in Fig. 6. Our results are as 
follows: 

— = 0.94±0.05, — = 0.98±0.02, — =0.92±0.01. 
This is compared with the result obtained using Wilson gauge action in Ref.[25]: 

— = 0.87 ±0.02 



Hence substantial improvement is obtained with the use of Iwasaki impoved action. 
This suggests the abelian monopole part could reproduce the full string tension in the 
continuum limit. 



Comparison of <B>/<A> Data 
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FIGURE 7. Non-abelianicity Q vs. Cooling Step, for several values of j3. This figure is from Ref.[10]. 



Abelian dominance in infrared QCD is also checked in some works with the use of 
controlled cooling. Giedt and Greensite [10] have measured the ratio B/A where 

B = -^I^-^IETrl^W,^)] 2 } 



>jm>n 



H P V X i>j 

and have shown that it decreases as the cooling steps as far as the string tension is 
kept non-vanishing. See Fig.7 taken from Ref.[10]. Cooling is expected to reduce the 
big Coulombic interaction in the non- Abelian case while keeping infrared confinement 
property. 

Cea and Cosmai[l 1] have measured connected correlation operators of Wilson loops 
and lxl plaquette using the controlled cooling and have obtained the penetration length 
in non- Abelian case. It is almost equal to the penetration length determined by Abelian 
Wilson loop and Abelian electric fields in MA gauge. See Fig. 8 taken from Ref.[l 1]. 

We also try to check Abelian dominance in Landau gauge using a controlled cooling 
[26] under which the string tension remains non-vanishing. For reader convenience let 
us, briefly, illustrate our cooling procedure. The lattice gauge configurations are cooled 
by replacing the matrix U^(s) associated to each link / = (s, p.) with a new matrix U'Js) 
in such a way that the local contribution to the lattice action 

S(s) = l-ltx{U fl (s)k(s)F(s)} (6) 

is minimized. F(s) = k(s)F(s) is the sum over the "U-staples" involving the link / and 

k(s) = y^det (P(s)j, so that F(s) E SU(2). In a "controlled" or "smooth" cooling step 
we have 

U^ S )^U^( S )=V(s)U^s), (7) 
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FIGURE 8. /I and /Ia (in units of A^-) versus j3 for square Wilson loops. Circles, squares, and triangle 
refer to L = 16, 20, 24 respectively. Crosses and diamond refer to the Abelian projected correlator p$, 
with L = 16, 20 respectively.This figure is from Ref. [11]. 



where V(s) is the SU(2) matrix which maximizes 

tr{V(s)U^s)F( S )} (8) 
subjected to the following constraint on the SU(2) distance between U^(s) and U'^(s): 
1 



-tr 



(ul{s)-U«{ S j) (UM-U^s)) 



(9) 



We adopt 8 =03. A complete cooling sweep consists in the replacement Eq. (7) at each 
lattice site. 

We find the profile of the non-Abelian electric field E% tends to that of the Abelian 
one of in the long-range region as shown in Fig. 9. This is consistent with the above 
result[ll] in a different approach. 



DIMENSION 2 GLUON CONDENSATE 

Now we have shown that the magnetic displacement currents are important in the dual 
Meissner effect when there are no monopoles. Then how can we understand the origin of 
the dual Meissner effect without monopoles? The Abelian dual Meissner effect indicates 
the massiveness of the Abelian electric field as an asymptotic field: 

{d 2 p -m 2 )E A ~0. (10) 

This leads us to a dual London equation which is a key to the dual Meissner effect. Let 
us evaluate the curl of the magnetic displacement current. Using Eq.(4), we get 



Vxd 4 B A = V(V.£ A )-V 2 £ A . 
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FIGURE 9. Profiles of non-Abelian electric field after cooling. stands for the Abelian electric 
field. W(R x T = 6 x 6) is used. 



From Eq.(lO), we get the dual London equation: 

V x d 4 B A ~ (d 4 2 - m 2 )E A . (11) 

Let us remember a simple mean-field approach developed by Fukuda[27]. Neglecting 
gauge-fixing and Fadeev-Popov terms, we have equations of motion D^F^ V = and 

the (non-Abelian) Bianchi identity D^*F b v = 0. Applying D operator to the Bianchi 
identity and using the Jacobi identity and the equations of motion, we get 

(D 2 p ) ab F b v = 2 g e abc F b a F$ a . (12) 

Notice {D 2 p ) ab = d 2 8 ab + g£ acb (d p A L p ) + g 2 (A a p A b p - 8 ab (A c p ) 2 ). Hence if < A a ^A b >= 
S ab 8^ v v 2 7^ 0, we see asymptotically that the electric fields become massive (d p — 
with m 2 = 8g 2 v 2 [28]. Now the Abelian electric field is also massive asymp- 
totically (dp — m 2 )E^ k ~ 0. Hence the dual London equation (1 1) is obtained. 

The importance of the dimension 2 gluon condensate has been stressed by Zakharov 
and his collaborators [13] and Refs.[14, 15]. Recent discussions on the value of the gluon 
condensate are seen in Ref.[16]. Some of them discuss the mass generation of the gluon 
propagator. But as pointed out by Fukuda[27], the non- vanishing dimension 2 gluon 
condensate leads us to a conclusion that field strengths F£ v (s) instead of gluon fields 
A^(s) become good canonical variables having a finite mass, whereas gluon propagators 

have a (p 2 )~ 2 behavior showing confinement. 

Although the operator of the gluon condensate is gauge-variant, it is proved recently 
the expectation value is gauge invariant[17]. Hence the gluon condensate has a physical 
importance, if the proof is correct. 
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